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ABSTRACT 

Recent work has shown that a young, rapidly rotating neutron star loses angular momentum 
to gravitational waves generated by unstable r-mode oscillations. We study the spin evolution 
of a young, magnetic neutron star including both the effects of gravitational radiation and 
magnetic braking (modeled as magnetic dipole radiation). Our phenomenological description 
of nonlinear r-modes is similar to, but distinct from, that of Owen et al. (1998) in that our 
treatment is consistent with the principle of adiabatic invariance in the limit when direct 
driving and damping of the mode are absent. We show that, while magnetic braking tends to 
increase the r-mode amplitude by spinning down the neutron star, it nevertheless reduces the 
efficiency of gravitational wave emission from the star. For B J> 10 14 (^ s /300 Hz) 2 G, where 
i/ s is the spin frequency, the spindown rate and the gravitational waveforms are significantly 
modified by the effect of magnetic braking. We also estimate the growth rate of the r-mode due 
to electromagnetic (fast magnetosonic) wave emission and due to Alfven wave emission in the 
neutron star magnetosphcrc. The Alfven wave driving of the r-mode becomes more important 
than the gravitational radiation driving when B 10 13 (^ s /150 Hz) 3 G; the electromagnetic wave 
driving of the r-mode is much weaker. Finally, we study the properties of local Rossby- Alfven 
waves inside the neutron star and show that the fractional change of the r-mode frequency due 
to the magnetic field is of order 0.5 (B/10 16 G) 2 (i/ S /100 Hz)~ 2 . 

Subject headings: gravitation - hydrodynamics - stars: magnetic fields - stars: neutron - stars: 
oscillation - stars: rotation 



1. Introduction 

Gravitational radiation drives unstable the r-modes of rotating neutron stars (Andersson 1998; 
Friedman & Morsink 1998) via the Chandrasekhar-Friedman-Schutz mechanism ("CFS instability"; 
Chandrasekhar 1970; Friedman & Schutz 1978a, b). Moreover, the r-mode instability, due to current 
moment radiation, is so strong that viscous dissipations present in young neutron stars are insufficient to 
suppress the instability (Lindblom, Owen, & Morsink 1998; Andersson, Kokkotas, & Schutz 1999; Lindblom 
& Mendell 2000). Owen et al. (1998) modeled the spin evolution of a newly-formed, rapidly rotating 
neutron star and found that the growth of the r-mode causes the neutron star to spin down by losing 
angular momentum to gravitational radiation and that the emitted gravitational waves can be a promising 
source for the LIGO/ VIRGO gravitational wave detectors currently under construction. 

It has been suggested that most neutron stars are born with relatively small rotation rates (much less 
than the maximum break-up rate) and moderate magnetic fields of 10 12 — 10 13 G (Spruit & Phinney 1998). 
These neutron stars are not efficient gravitational wave emitters since the r-mode instability is suppressed 
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by viscosity at slow rotations. However, recent observations of soft gamma-ray repeaters and the so-called 
"anomalous x-ray pulsars" (e.g., Vasisht & Gotthelf 1997; Kouveliotou et al. 1998, 1999) indicate that there 
exists a distinct population of neutron stars with much stronger magnetic fields, i.e., magnetars (Thompson 
& Duncan 1993, 1995, 1996). The critical factor that determines whether a proto-neutron star becomes an 
ordinary pulsar or a magnetar may be its initial spin. It is conceivable that magnetars are born with rapid 
rotations which facilitate dynamo actions that generate their superstrong magnetic fields. It is therefore 
possible that the r-mode instability may only be relevant to proto-magnetars. 

In this paper we examine several issues related to the magnetic field effects on the r-modc and on 
the spin evolution and gravitational wave emission of a newly-formed, rapidly rotating neutron star. An 
obvious issue is the magnetic braking of the neutron star due to magnetic dipole radiation. While magnetic 
braking is negligible compared to r-mode-driven spindown for canonical magnetic field strength (~ 10 12 G) 
and rapid rotation, it becomes important for B <; 10 14 G and slower rotation. Magnetic braking can also 
drive/maintain the r-modc instability indirectly by changing the background state of the star according to 
the principle of adiabatic invariance. These effects of magnetic braking on the spin and r-mode evolution, as 
well as the gravitational waveforms, are studied in §2. In §3, we discuss several magnetic effects on r-modes, 
including r-mode instabilities driven by electromagnetic (fast magnetosonic) wave emission and by Alfven 
wave emission in the neutron star magnctosphere. Section 3 also includes a discussion of how magnetic 
fields modify the r-mode frequency and structure. We summarize our results and discuss their implications 
in §4. 



2. Spin and R-mode Evolution 
2.1. Basic Equations 

We first summarize the basic properties of r-modes needed for our discussion [see Owen et al. (1998) 
and references therein] . To leading order in Q s (the angular spin frequency of the neutron star) , the Eulerian 
velocity perturbation of the r-mode is of the form Sv = aQ s R(r / R^Y^ e lwt , where a is the dimensionless 
amplitude of the mode, R is the neutron star radius, = [1(1 + l)] _1 / 2 rV x (rVY; m ) is the magnetic-type 
vector spherical harmonic of order l,m. Throughout this paper, wc consider only the dominant I = m, = 2 
mode. The angular frequency of the mode in the inertial frame and the absolute value of the corresponding 
frequency in Hz are w = — 4f2 s /3 and v = |w|/(27r) = 4z^ s /3, where v s is the spin frequency (Papaloizou 
& Pringle 1978; Provost, Berthomieu, & Rocca 1981; Saio 1982). The frequency of the gravitational wave 
emitted by the r-modc is also v. We shall adopt the same neutron star model as in Owen et al. (1998), 
i.e., M = 1.4M Q , R = 12.53 km, and polytropic index r = 2. For such a neutron star, the damping 
rate of the mode (i.e., in the linear regime, the mode amplitude evolves as e ^*-*/T) can 

be written as 

t -1 = tgr -1 + TV -1 , where tqr is the gravitational radiation timescale as given by 



1 ^fA^GMR^ I 
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and Ty is the viscous timescale as given by 
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Here Cl s = n s /y^Gp, f GR = -3.26 s, f s = 2.52 x 10 8 s, f B = 6.99 x 10 8 s, J = 1.635 x 1CT 2 , and 
T = 10 9 T g K and p are the temperature and mean density of the neutron star. 
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In the presence of r-mode oscillations, the neutron star will be differentially rotating. As gravitational 
radiation takes away angular momentum from the r-mode, the differential rotation is expected to grow.[] 
Viscous forces tend to erase the differential rotation, i.e., to suppress the mode growth, by transferring 
angular momentum between the mode and the mean stellar rotation. In the absence of a precise nonlinear 
description of r-modes, we shall derive the equations governing the evolution of the rotation frequency and 
the r-mode amplitude using the phenomenological approach similar that of Owen et al. (1998), although 
our equations differ from theirs in an important aspect, even when B = 0, as we explain below. 

The total angular momentum of the neutron star can be decomposed into the angular momentum due 
to a mean rigid rotation, J = Itt s , where / is the moment of inertia, and the canonical angular momentum 
of the mode,F| 

3 ~ 

J c = --JMR 2 n s a 2 . (2.3) 

The canonical angular momentum of the mode can increase through gravitational radiation and decrease 
by transferring angular momentum to the star through viscosity. Thus 

dJ c 2 2 , , 

-J- = Jc Jc + Mother, 2.4 

at t gr t v 

where Aether is the torque associated with other processes that can drive the mode. Section 3 discusses 
two possible mechanisms that can contribute to -/V thcr, namely the driving of the mode by electromagnetic 
radiation and by Alfven wave emission. In this section, we consider driving only by gravitational radiation, 

SO that Aether = 0. 

Another evolution equation is obtained from noting that the angular momentum associated with 
the mean rigid rotation can increase through a transfer from the mode due to viscosity and decrease by 
magnetic braking. Thus 

f^ = -l Jc _^, (2 . 5 ) 
at Tv tm 

where {—I$I s /tm) is the magnetic braking torque and tm is the corresponding timescale. Adopting the 
simplest magnetic dipole model, we have 

1 _ B 2 R 6 Q 2 _ 1 2 A2 _ 1 2 ( \ 2 

m~ 6cH ~ ^ Blis ~ 1.21 x Kp s Bli \Tttb) • (2 - 6) 
Here tm = 6.74 x 10 4 s, B = 10 14 £?i4 G is the neutron star magnetic field, and we have used / = IMR 2 



with / = 0.261 for T = 2 polytrope. Substituting equations (2.3) and (2.6) into equation ( p.q ) yields 



— = -2Q^ ±, (2.7) 

at tv tm 



1 This is borne out by the exact solution of the evolution of an unstable bar-mode in an incompressible Maclaurin spheroid 
as it evolves toward a Dedekind ellipsoid driven by gravitational radiation reaction (Miller 1974; Detweiler & Lindblom 1977; 
Lai & Shapiro 1995). See also Spruit (1999) and Levin & Ushomirsky (2000). 

2 The physical angular momentum associated with the r-mode, denned as the change in angular momentum of the star due 
to the mode, generally differs from J c (Friedman &; Schutz 1978a; Levin & Ushomirsky 2000). This difference amounts to a 
proper definition of the mean rotation rate Q B for the differentially rotating star. The usefulness of J c lies in the fact that in 
the absence of radiation and viscous dissipation, J c is conserved on account of the first-order perturbation equations, and that 
J c = —m l c is proportional to the wave action I c , an adiabatic invariant (Friedman & Schutz 1978a); see discussion following 



equation (2.8) 
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where Q — 3J/2I = 9.40 x 10 2 . From equations (2.4) and (2.7) we obtain the evolution equation for the 
mode amplitude, 



da 



TGR 



TV 



2r M 



(2.8) 



Note that our evolution equations (2.7) and (|2.8|) differ from those given in Owen et al. (1998) even for 



B = 0. The difference arises because Owen et al. (1998) derive one of their evolution equations by analyzing 
the rate of change of the mode energy (in the rotating frame) , namely 



dE c 
dt 



-IE, 



(- 



\TGR Ty 



(2.9) 



where E c = JMR Qjcf /2. This is in contrast to our equation (2.4) (with N ( 



'other 



0) which is based 



on angular momentum considerations. In the regime where the spin frequency does not change, the two 
treatments give the same result for the evolution of the mode amplitude, since E c = (— cr/2) J c , where 
a = 2Q s /3 is the mode frequency in the rotating frame. However, when £! s changes, the two treatments 
yield different results. We believe that equation ( [2.9| ) is conceptually incorrect in the regime when both 
Q s and a are varying. This can be most easily understood by considering a hypothetical situation where 
there is no gravitational radiation and viscosity, but an external torque, such as magnetic braking, slows 
the star down gradually. In such a situation, an initial wave, such as a r-mode, that propagates on the star 

— 1/2 

will be amplified in amplitude by a oc fl s according to the conservation of wave action (an adiabatic 



invariant), which is simply proportional to J c . On the other hand, equation (2.9) would indicate a oc Q 
for r, 



GR 



0, which is clearly incorrect. Nevertheless, the difference between our evolution equations 



and those of Owen et al. (1998) comes in at the level of a 2 Q ~ 0.1a 2 , which is small unless the mode 
becomes strongly nonlinear (a > 1). However, when the magnetic braking effect is comparable to or 



greater than the gravitational radiation effect, adopting equation (2^9) would give a significant error in the 
predicted behavior of the mode amplitude. 

There are three different stages in the evolution of the r-mode. In the first stage, when the mode is 



still in the linear regime (with a <J 1), the mode amplitude grows according to equation (2.8), while the 



spin evolves according to equation (|2.7|) . Note that for <C |tqr| , the mode grows exponentially in 
this stage, but for t^ 1 ^ |tgr| _1 , the mode growth can be significantly slower (see below). The second 
stage occurs when nonlinear effects prevent further growth of the mode. We adopt the same ansatz as in 
Owen et al. (1998), i.e., when the mode amplitude reaches a sa t = yfti ~ 1, it will be held constant, so that 
da/dt = 0. In this stage, the spin evolution is then described by 



dCl s 
~dt 



2a 



i 



t gr 1 - kQ t m 1 - kQ 



(2-10) 



Finally, the third stage occurs when the star has been slowed down, and viscosity starts to damp out the 
mode. This stage begins when the right-hand side of equation (2M) becomes negative. In the third stage, 



the evolution equations are again given by equations (2.7) and (2 



To evaluate the viscosities, we adopt the same temperature evolution of the neutron star, based on 
neutrino cooling via the modified URCA process, as in Owen et al. (1998), i.e., T = 10 9 K(i/l yrT 1/6 for 
t ;> 10 s. [For t <, 10 s, T ~ 10 10 — 10 11 K, the star is opaque to neutrinos, and this bulk viscosity expression 
is invalid.! 
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2.2. Results 



Integrating equations ( |2.7D , (2.8), and ( p.lOj ) numerically, we can easily obtain the mode amplitude 



and neutron star spin as functions of time. Figure [l] gives examples of such evolution for several different 
magnetic field strengths: B = G, 10 14 G, 10 15 G, 1.2 x 10 15 G, and 2 x 10 15 G. We begin the star 
rotating at Ct S Q — &k — 2/3, which corresponds to v s q = 890 Hz. We choose the initial mode amplitude 
«o = 10~ 4 and the saturation amplitude a sa t = 1. In the absence of magnetic braking (B — 0), the mode 
grows exponentially during the first stage and reaches saturation at t sa t — 350 s, while the change in the 
spin frequency during this stage is on the order of 0.1 Hz (Owen et al. 1998); the neutron star then spins 
down after 1 year to v s ~ 90 Hz. We see that for relatively weak magnetic fields (B <, 10 14 G), magnetic 
braking affects the spin evolution only at later times, when the mode has saturated and when Tj^ 1 cx Bf 4 Q 2 
becomes comparable to |tgr| _1 oc fl®. However, for sufficiently high magnetic fields (B <; 10 15 G), the 
growth of the mode can be significantly delayed (e.g., t sat ~ 1700 s for B = 10 15 G), and the star spins 
down appreciably even in the first (linear) stage. 

It is of interest to note that in the presence of magnetic braking, the mode does not begin damping out 
even after the stellar spin has dropped to well below 100 Hz, the threshold spin frequency for the r-mode 
instability without magnetic braking. This comes about because magnetic braking can drive (albeit slowly) 
and maintain the oscillation mode by changing the background state, i.e., rotation, of the star in which the 



mode lives, in accordance with the principle of adiabatic invariance. Indeed, equation (2^8) indicates that 
the critical spin frequency for r-mode growth is determined by the condition 

1 -U^>0, (2.11) 



|tgr| tv 2 r M 

where we have set (l — a 2 Q) ~ 1. Figured depicts the critical spin frequencies as functions of the neutron 
star temperature for several different values of B. The evolutionary tracks are also shown for comparison. 
It is evident that magnetic fields decrease the requisite critical spin frequencies for mode growth by 
countering the viscous damping. At late times, when the neutron star has cooled, it is the shear viscosity 
that damps the mode. Without magnetic braking, this viscosity prevents the growth of the mode when 
v s < 65 T g Hz. However, magnetic braking can sustain the growth of the mode if v a > 31 T^ 1 B^l Hz. 

The effects of magnetic braking depend on both B and the initial spin f2 s o of the neutron star. The 
critical relevant parameter is the ratio of the magnetic braking timescale and the gravitational radiation 
timescale, 

P=^ = — (-X^flA = 0-645- f • (2.12) 

|t G r| 75 \3J c 4 B 2 14 Vl00Hz/ v ' 

If the initial spin Q s q is such that |tgr| _1 Ty 1 or |tm| _1 S> Ty 1 , we can neglect the viscous term in 



equations (2.7) and fl2.8p. We then find that i sat and r2 S at, the time and rotation rate at the saturation 



point, respectively, are determined by 

1 = -^;lna; sat + -^(1 -a; 4 at ), y aat = i sat /t sat o = ^(l/a&t _ !)j ( 2 - 13 ) 

where /3 is the initial /?, x sat = r2 sa t/O s o, 7 = In (a sat /o;o), and i sat0 = |t G r| In (a sat /ao) is the time to 
saturation without magnetic braking. For P0/7 3> 1, we find 

^at = l-^-£f|-^-V 1^ = 1+^3-^^. (2.14) 



27/ Po 
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Figure || shows the general results for a; sa t and y sat with 7 = lnlO 4 and lnlO 8 . It is clear that magnetic 
braking increases both i sat and Af2 sat = il s o — f2 Ba t = (1 — ^sat) ^so- Note that in the B = limit, Af2 sat 
is determined by viscosity and is of order 10 -4 ; the expressions above are valid as long as the resulting 
AO sa t is greater than the zero-field limit. From Fig. ||, we see that the linear evolution will be significantly 
affected when 

A) £ 5 7 ~ 50 v s0 £ 300 (^) B\{ 2 Hz. (2.15) 

The transition is quite dramatic. For example, we see from Fig. [l] that the time it takes for the mode to 
saturate is t sat ~ 1700 s for B = 10 15 when v s0 = 890 Hz. But for B = 1.2 X 10 15 G, saturation has not 
been reached even after t = 1 yr. This will significantly affect the gravitational wave signals (see §2.3). In 
the second (saturation) stage, the spindown rate is given by [see eq. ( 2. 10] )] 



]_d^s k f v s Y (1 1 \ (216) 

v a dt 3yrV100Hz/ \ 2nQ(3 J ' { ' ' 

Thus, even for (3q ^> 50, so that the linear evolution is not affected, magnetic braking can still affect the 
evolution in the the second (saturation) stage, when 

/3£-L~5 v s < 170 k^bH 2 Kz. (2.17) 

ZkLJ 



2.3. Gravitational Waveforms 

The rate of angular momentum loss due to the r-mode can be written as 

dJ c \ 2 T lu 2 , 9 , . . 

-jf) = J c = 4nd 2 -—(h 2 + + h 2 x ), (2.18) 

dt J GR t G r uj 167T 

where h+ and h x are the strain amplitudes for the two polarizations of the gravitational wave and d is the 
distance to the source (Shapiro & Teukolsky 1983). Following Owen et al. (1998), we define the averaged 
gravitational wave amplitude h(t) by [h(t)] = (3/10) (h\ + h 2 x ). Using equations ( |2.3| ) for J c and (2.1) for 
tgr i we find 

Kt) = 2 -§M^^ - 1-8 x 10"-a (-^-) 3 ( ^ . (2.19) 
45 V 30 c 5 d VI kHz/ \ d J y ' 

The bottom panel of Fig. [I] plots h(t) as a function of time for several different magnetic field strengths, 
with the distance to the source d — 20 Mpc. As expected, the effect of magnetic braking is to decrease the 
gravitational wave amplitude and shorten the duration of gravitational wave emission. 

To examine the detectability of gravitational waves, we define the characteristic gravitational wave 
amplitude h c via 

h ^^^ = v WJwr- (2 - 20) 



where h is the Fourier transform of the gravitational wave amplitude. In equation ( 2.20 ), dv/dt can be 



obtained from equations (2.7) and (2.10) for the three stages of spindown evolution. Figure [I] shows h c {v) 



for two different initial spin frequencies (y s o = 890 Hz and 400 Hz), each with several different values of B. 
The other parameters are the same as in Fig. [j] and Fig. ^. For relatively weak magnetic fields, the vertical 
spikes seen in Fig. [| result from the fact that when /?□ S> 50, the spin frequency, and hence v, is essentially 
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constant before saturation.^ For higher magnetic field (smaller the spikes are broadened, and the width 
of the spike, Av = (2/3tt)(1 — x sa _t)^ s o, becomes significant for /? ^ 50 [see eq. ( 2.15 ) and Fig. Near 
saturation, the gravitational wave amplitude is 

h c {u = , sat ) = 4.9 x ID"" (j^-y 1 "' \ 1/2 *llhl£ fl + 118^ V2 ■ (2.21) 



For small B, we have x sat = y sa t = 1, and equation (2.21) reduces to the zero-field result. In the saturated 
regime, we have 

M ,).5.8xl0--(^) 1/2 (l + ^) . (2.22) 



Thus for [3 <y 1/(2kQ) ~ 5 [see eq. (2.17)], the characteristic amplitude is reduced from the zero-field value. 



In Fig. [|, we also compare h c (u) with the expected RMS noise spectra h lms = \J vSyJy) of three (initial, 
enhanced, and advanced) generations of LIGO, where Sh (v) is the power spectral density of the detector 
strain noise. The expressions for Sh(v) are given in Owen et al. (1998). We see that the signal-to-noise 
ratio, proportional to v~ 1 {h c /h rras ) 2 integrated over the frequency band, can be significantly reduced by 
the magnetic braking effect. 



3. Effects of Magnetic Fields on R-modes 

Magnetic braking discussed in §2 only indirectly affects the r-mode instability. In this section, we 
discuss several direct effects of magnetic fields on r-mode oscillations and their instabilities. 



3.1. Driving of R-Mode by Electromagnetic Radiation 

R-mode oscillations in the neutron star perturb the magnetic field anchored on the star, and this 
leads to electromagnetic radiation. Just like gravitational radiation, the electromagnetic radiation also 
drives r-modes unstable via the CFS mechanism. Consider the I = m = 2 r-mode oscillation with surface 
Lagrangian displacement 

^olV^lr^We-^ (3.1) 



sin0<9<£' 86, 

where u) = 40 s /3 (we shall adopt the convention lu > in this section and in Appendix A). The power of 
electromagnetic radiation associated with this oscillation is given by (see Appendix A) 



EM 



cR 2 ( w£ n \ f 



1087T V c J \ c 



B — . (3.2) 



Since the mode energy (in the inertial frame) is Ec = ~(9/56tt)Mu> 2 ^q for a uniform star, the timescale 
of EM driving is then 

' 7 /1J2fl " ' ' //V ^ = U r Afu^. (3.3) 



t em 3 5 V Mc 2 J \ c j (3.8 x 10 6 s) 



3 Recent work by Levin &; Ushomirsky (2000) suggests that the spindown evolution may be unrelated to mode saturation. The 
mode frequency in this case may not be constant before saturation, and thus the vertical spikes in Fig. ^ would be broadened. 
The details, however, are uncertain at this time. 
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Comparing this with t g ^ = — (4 8 /3 7 )(GMi? 4 i7g/350c 7 ) (for a uniform neutron star), we find 

p™^A7M^Rl BlJ^)\ (3.4) 

1/TGR \ V s J 

where M\a is the neutron star mass in units of 1.4Mq, R\q is the radius in units of 10 km, and B\q is the 
magnetic field in units of 10 16 G. Clearly, electromagnetic driving of r-mode is negligible unless the neutron 
star has B ^ 10 16 G. 

Since a real neutron star possesses a magnetosphere, whose plasma frequency is much greater than 
the typical rotation frequency, it is unlikely that electromagnetic radiation can be emitted from the star. 
However, as pointed by Blaes et al. (1989), in a relativistic plasma, fast magnetosonic waves have the 
same dispersion relation as electromagnetic waves, and their emission is therefore analogous to that of 
electromagnetic waves. 



3.2. Driving of R-Mode by Alfven Wave Emission 

The shaking of magnetic field lines by r-mode oscillations also results in Alfven wave emission to the 
magnetosphere. The torsional Alfven waves carry away energy and angular momentum from the r-mode and 
thus drive it unstable. Because of the uncertainties in the magnetosphere physics, here we shall be content 
with an order-of-magnitude estimate of the Alfven wave power. In a relativistic plasma, the group velocity 
of Alfven waves is close to the speed of light and is along the magnetic field line. A fluid displacement £o 
on the neutron star surface produces a kink in the attached field line which in turn attempts to straighten 
up by launching Alfven waves. In one oscillation cycle, the signal travels a distance ~ c/oj. Therefore the 
amplitude of the Alfven wave is of order SB ~ (uj£o/c)B. Note that this is smaller than the magnetic 
field perturbation, of order (£o/R)B, inside the stellar surface. The power radiated by the Alfven waves is 
then Pa ~ [(SB) 2 / (S-k)](At:R 2 )c = (Rlo^qB) 2 /(2c). An improved calculation, taking into account various 
geometric factors, yields 



i ! /? ) \ B 2 ^ 2nR 2 \B-f \ sin 6>d6>~ 0.04 cR 2 \ —B ) . > 3.:. ) 



where we have assumed the field to be dipolar outside the stellar surface.^] Obviously equation ( |3.5|) should 
be considered as an estimate since we do not solve for the global field perturbation self-consistently. The 
timescale of Alfven wave driving is then 

^~-°- 4 (^)(5)=-2^ M ^» B »- ("» 



4 Note that the power rad iate d by Alfven waves scales as lu 2 while the power emitted by electromagnetic (fast magnetosonic) 
waves scales as io 6 [see eq. ( h.2| . The weak coupling of the magnetosonic waves can be understood in the following way [see 
Blaes et al. (1989) for a similar discussion on elastic waves]. Prior to reaching the stellar surface, the oscillation mode has 
the dispersion relation w ~ kv, where v ~ wR. The magnetosonic and Alfven waves in the magnetosphere have the dispersion 
relation w = kc and u> = kecostp, respectively, where rjj is the angle between k and B. From the continuity of the projection of 
the wave vectors onto the incident boundary, ki sin 8 j = ktsinBt (Snell's law), where 9i and 9t are the angles of incidence and 
transmission, we have s'mOt = (c/v) s'm8i and sinSt = (c/v) cosi/>sin0i for the magnetosonic and Alfven waves, respectively. 
Since c>«, the oscillation modes must have a very small angle of incidence in order for them to be transmitted (coupled) to 
the magnetosonic waves, otherwise the mode is reflected at the boundary. On the other hand, the additional degree of freedom 
allowed by cos ip permits a wide range for the angle of incidence of the oscillation modes in the case of Alfven waves, so that 
the modes are more easily coupled to Alfven waves, with the only restriction that k be almost orthogonal to B. 
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Comparing the Alfven wave power to the gravitational radiation power, we have 



1/ta « 999 /100Hz\ 6 , s 

' 10 3 M^Rlo 2 BfJ — ) . (3.7) 



l/ r GR V 

1 /3 1 /3 

Thus, when the spin frequency v s <> 320 Hz B-^ = 150 Hz £^3 , Alfven wave emission becomes more 
important than gravitational radiation in driving the r-mode instability. 



We note that our calculation for the Alfven wave power, equation ( p.5[ ), includes the energy transmitted 
through both the open and closed magnetic field lines. The energy transmitted through the open field lines 
will be dissipated at infinity; the energy associated with the closed field lines is less clear. Also, only the 
coherent part of the Alfven wave carries away canonical mode energy, thereby driving the mode. In the 
absence of a complete understanding of the global properties of the magnetosphere, we cannot definitively 
calculate the mode driving/damping rate due to Alfven wave emission. If we include only the power radiatied 



through the open field lines, then equations (3.5)-(3.7) must be modified. The integral in equation is 



instead evaluated from 9 = to 9 = 9 cap = sin 1 ^Q s R/c, This results in the power and driving rate of 
Alfven wave radiation being multiplied by a factor of ~ 4 sin 4 9 cap ~ 4(fi s i?/c) 2 ~ 2 x I0~ 3 i? 2 (z4/100 Hz) 2 . 

An important caveat in our calculations of ta and tem is the neglect of neutron star crust, which is 
expected to form for T <, 10 10 K. The r-mode oscillation in the fluid core may not penetrate the crust. As a 
result, the r-mode-induced perturbations to the magnetosphere are greatly reduced except for hot neutron 
stars. 



3.3. Rossby- Alfven Waves 

A r-mode oscillation perturbs the magnetic field inside the star, and the resulting magnetic stress acts 
as a restoring force for the oscillation. If the magnetic field is sufficiently strong so that the Alfven speed 
becomes comparable to the phase velocity of the r-mode, the oscillation frequency and eigenfunction will 
be significantly modified by the magnetic field. In essence, the oscillation does not correspond to a pure 
Rossby wave but contains a mixture of Rossby wave and Alfven wave characteristics. 

While a global analysis of the Rossby- Alfven mode of a rotating magnetic star is beyond the scope 
of this paper, it is instructive to examine the properties of local Rossby-Alfven waves. We consider fluid 
perturbation £ in the neighborhood of (r, 9, (f>) inside a fluid star which rotates at angular frequency 
Q s around the z-axis. It is convenient to set up a local Cartesian coordinate [x\, x%, X3) so that 
dx\ = rd9, dx2 — r sin. 9d<f>, dx% — dr. We consider the simplest case in which the unperturbed magnetic 
field B is along the radial direction, and £ has only a; 1,2:2 components. The oscillation involves no density 
perturbation and thus satisfies V • £ = 0. The perturbed magnetic field is given by 

SB = V x x B) = B^-l (3.8) 

where we have neglected the curvature of the unperturbed magnetic field. The x%, X2 components of Eider's 
equation in the rotating frame can be written as 

^ = JJJH + 2(l , a 4 + f^ (3.9) 
at 2 p dx\ ot Airp dx^ 



<9 2 6 1 dSP ^ 06 t B 2 2 £ 

4 



dt 2 ~ p dx 2 ^ dt + 4vrp fe 2 ' (3 ' 10) 
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where tt r = f2 s cos 9, and the last terms of the equations represent the magnetic force per unit mass given 
by Fb = (B ■ V)(5B/(47r j o). (The displacement current in the Maxwell equation is neglected since the Alfven 
speed inside the star is much less than the speed of light.) We can eliminate the pressure perturbation 6P 
from equations (3J3) and ( |3.10| ), and then set £ oc e l,yt + lk i x i+' lk 2X2+ik r x 3 ^ wnere a j s ^he mo de frequency 



in the rotating frame. Combining the resulting equation with V • £ = + fe^) = 0, we obtain the 

dispersion relation 



where ua = (B 2 /Anp) 1 / 2 is the Alfven speed. For B = 0, this is simply the dispersion relation for local 
Rossby waves. Indeed, with dfl r /dxi = — O s sm0/r, &2 = m/(r sin ff) and k\ + fcf = 1(1 + l)/r 2 , the 



right-hand-side of equation ( 3.11 ) becomes Co = 2mCl/[l(l + 1)], the r-mode frequency of order l,m. For 



B ^ 0, the Rossby- Alfven wave frequency is given by 

a = ^±i(a 2 +4 CT |) 1/2 , (3.12) 
where cta = WAfc r is the local Alfven angular frequency 

ct a ^ 300 (k r R)B 16 p~4 /2 R^ s _1 . (3.13) 
For ctq 3* CAj the change of the r-mode frequency, with the I = m — 2 mode and k r R ~ 1, is of order 

^ = ^0.5 M r„ ! (^V. (3.14, 

For Aa/ao <C 1, the oscillation manifests itself as Rossby waves in the bulk interior of the star, and the 
magnetic energy associated with the oscillation is negligible. But near the stellar surface, where a a ^> er , 
the oscillation assumes the Alfven wave character. 



4. Discussion 

We have studied the combined effects of magnetic braking and gravitational radiation generated 
by unstable r-mode oscillations on the spin evolution of a young, rapidly rotating neutron star. Since 
the gravitational radiation timescale of r-mode scales as ^ s -6 , where v s is the spin frequency, while the 
magnetic braking time scales as (Bv s )~ 2 , it is obvious that the magnetic braking becomes important only 
for sufficiently large B and small v s [see eqs. (2.15) and (2.17)]. We have examined the characteristics of the 
gravitational waveforms for different magnetic field strengths. As discussed in §1, it is possible that young 
rapidly rotating neutron stars are necessarily endowed with superstrong magnetic fields (B ^> 10 14 G); thus 
one (an optimist) might suggest that the spectral features of the gravitational waves can be used to probe 
the initial magnetic fields and spins of these neutron stars. 

Magnetic fields can affect the r-mode oscillation in different ways. We have discussed a few of these 
effects in this paper (see §3). For sufficiently large B and slow rotation, the driving mechanisms of the 
r-mode due to electromagnetic (fast magnetosonic) wave emission and Alfven wave emission can compete 
with the driving due to gravitational radiation. In an interesting paper, Rezzolla et al. (2000) suggested 
that r-mode oscillations can wind up the stellar magnetic field due to differential drifts of fluid elements^] 



5 We note that in the case of an incompressible Maclaurin spheroid, where the nonlinear m = 2 f-mode is exactly known and 
is described by a Riemann S-ellipsoid, there is no differential drift and therefore no secular winding of field lines. 
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and this serves as a damping mechanism for the r-mode. The damping rate estimated by Rezzolla et 
al. (2000), if correct, is about (c/Cl s R) times larger than our estimated growth rate due to Alfven wave 
emission (assuming that the toroidal magnetic field is of the same order as the poloidal field). It is of 
interest to note that their estimate of the damping rate, 7, is of the same order of magnitude as the change 
in the mode frequency due to the magnetic field, i.e., 7 ~ Act ~ g\I<j§ (see §3.3). 

Bildsten & Ushomirsky (2000) have shown that the presence of a solid crust in the neutron star can 
strongly enhance the viscous damping of the r-mode.[] For viscous damping 10 5 times stronger [as suggested 
by Bildsten & Ushomirsky (2000)1, the critical spin frequency below which the mode begins to damp is 

— 1/3 

given by v s ~ 440 T g 1 Hz for low magnetic fields. It is of interest to note that when the magnetic field is 
sufficiently strong, the mode does not decay until v s < 100T 9 _1 i3 1 ^ 1 Hz due to the indirect driving of the 
mode by magnetic braking [see eq. ( f2.ll )]. Nevertheless, for B <; 10 15 G, the strength of the gravitational 



waves is suppressed because magnetic braking dominates over gravitational radiation in the spindown, as 
shown in Fig. ||. 



Some of the calculations presented in this paper were done in the summer of 1998 while D.L. was at 
the Aspen Center for Physics. He thanks the participants of JENAM'99 "Gravitational Wave Astronomy" 
Workshop (September 1999, Toulouse, France) for their interest in this work, which prompted us to write 
this paper. He also acknowledges helpful discussion with Ira Wasserman, who has independently considered 
some of the issues addressed in this paper. Finally, we thank Yuri Levin for useful comments. This work is 
supported by NASA Grant NAG 5-8484 and an Alfred P. Sloan Foundation fellowship to D.L. 



A. R-Mode Driving by Electromagnetic Radiation 

Consider a star with a uniform magnetic field. We set up a coordinate system (x,y, z) centered on the 
star with the z-axis along J7 S , such that fi s = fl s z. Let the angle between the z-axis and the magnetic field 
axis be a. The magnetic field is then given in spherical coordinates as 

B = Bq [cos a cos 6* + sin a sin0 cos(0 — fi s i)] f 

+Bo [— cos a sin 9 + sin a cos 8 cos (cf> — Q s t)] 8 

-Bq [sina sin(<^ - Q s t)] (Al) 



The Lagrangian displacement of the r-mode at the stellar surface is given by equation (3.1), so that the 
velocity perturbation is 

"--^(s^-tt')'"**- ,A2 » 



6 The boundary layer analysis of Bildsten & Ushomirsky (2000) hinges upon the existence of a sharp (width less than a 
few cm) transition between the liquid core and the inner crust, i.e., a first order phase transition. While this may well be 
true (Baym et al. 1971), one cannot be completely sure. Note that the pressure scale height of this region is of order 10 4 cm. 
It is also not clear how rigid the the core-crust transition region is, given that the strong interaction plays a more dominant 
role than the electrostatic interaction in determining the properties of the medium, contrary to the situation in a classical 
one-component plasma for which the original melting criterion is obtained. Nevertheless, it clear that the crust can enhance 
the viscous dissipation, although one might consider the estimate of Bildsten & Ushomirsky (2000) as an upper limit. Recent 
calculations by Andersson et al. (2000) suggest that Bildsten & Ushomirsky (2000) overestimated the strength of the damping. 
Owen (1999) and Andersson et al. (2000) also point out that the energy released by a mode that is already present may delay 
the formation of the crust, thus negate the large viscous dissipation due to the boundary layer, and allow the spindown evolution 
to proceed by the method described in our paper. 
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Note that, for the discussion here, we have changed the sign of uj so that uj = 4f2 s /3 > 0. The star then has 
an induced electric field E = (<5v/e) x B, which possesses three components with different frequencies. 
These are 

E roodeil = -Mucosa cosd(^§+^Y 22 $)e-^ (A3) 



E -*|^» sina e -*-*(<-fc)« ,A4, 



E modc , 3 = -*^ 8 in (w^ + ^) e - , H^)*. (A5) 

Thus the I = m = 2 mode gives rise to electromagnetic waves with frequencies uj, uj + Q s , and uj — f2 s . Using 
the method and notation described of Jackson (1975), we can obtain coefficients a~E,(l, m) and ay[(l, m) which 
specify the electric {l,m) multipoles and magnetic (l,m) multipoles. For the electromagnetic radiation 
corresponding to ( |A3| ), the nonzero coefficients, up to I = 3, are 

, 00 , 2 (kR) 3 iutoBp 

a M (22) = -= cos a (A6) 

V6 3i c 

1 / - ( 1 \ (kR) 5 iujCoBo 



a E (32) = ^+ j__ —cosa, (A7) 



where k = uj/c. For the electromagnetic radiation corresponding to (A4), 



(A8) 



f11l 1 ( Rj. 1 ^ (^) 3 !^Mo . 
a E (ll) = -= V5H F ; smo 

x 2 (kR) 3 iujfnBn , , . „. 

a M 21 = ^sina A9 

V6 3i 2c 

flE(31) = 7!1V35^ 2c- Sma ' (M0) 



where k ~ (uj — Q s ) jc. For the electromagnetic radiation corresponding to (A5), 



,„„x 8 h(kRfiuj^Bo . , A11 , 
«e(33) = -^-^ 2^ Sma ' (AU) 



where fc = (w + f2 s ) /c. To leading order in (ujR/c), the contribution to the power from (A3) and (A4) are, 
respectively, 

Pi = ^^^^(cosa) 2 (A12) 
fil 

p > = s^^ B » <si " Q) • < A13) 

The power from ( |A5| ) is negligible. Comparing the total power emitted from electromagnetic radiation to 
the rate of energy loss due to gravitational radiation gives the timescale for electromagnetic driving of the 
mode in terms of the timescale for gravitational driving, namely 

^ = 0.47 M^RlA (cos 2 a + | sin 2 a) (^) ' . (A14) 
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Fig. 1. — Spin frequency v s , mode amplitude a, and gravitational wave amplitude h as functions of time for 
a 1.4 Mq, 12.53 km neutron star at d = 20 Mpc with v s q — 890 Hz, ao = 10~ 4 , and k = 1. The dashed 
line is for the spindown caused by just gravitational radiation (.6 = 0). The solid lines are for the spindown 
that also includes magnetic braking with B = 10 14 G, 10 15 G, 1.2 x 10 15 G, and 2 x 10 15 G from top to 
bottom, respectively. Note that the line for h with B = 2 x 10 15 G is not shown because it is below the level 
of interest. 
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Fig. 2. — Spin frequency v s as a function of the temperature T for a 1.4 M Q , 12.53 km neutron star 
with v s q = 890 Hz, ao = 10~ 4 , and k = 1. The thick long dashed line is for the spindown caused by 
just gravitational radiation (B = 0). The thick short dashed line and thick solid line are for the spindown 
that also includes magnetic braking with B = 10 14 G and B = 10 15 G, respectively. Also shown are the 
critical spin frequencies above which the r-modes are unstable to growth. These are denoted by thin lines 
corresponding to the different magnetic field strengths. 
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Fig. 3. — x S at and y sa t as functions of /?o/7 for a 1.4 M , 12.53 km neutron star with u s q = 890 Hz. The 
solid line is for a = 10~ 4 , and the dashed line is for a = 1CP 8 . See §2.2. 
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Fig. 4. — Characteristic gravitational wave amplitude h c — vh as a function of frequency f for a 1.4 Mq, 
12.53 km neutron star at d = 20 Mpc with v s0 — 890 Hz and ^ s0 = 400 Hz, a = 10 -4 , and n = 1. The 
short dashed line is for the spindown caused by just gravitational radiation (.6 = 0). The solid lines are for 
the spindown that also includes magnetic braking. For ^ s o = 890 Hz, the magnetic fields are 6 = 10 14 G, 
4 x 10 14 G, 8 x 10 14 G, 1 x 10 15 G, 1.1 x 10 15 G, and 1.2 x 10 15 G from top to bottom, respectively. For 
z/ s0 = 400 Hz, the magnetic fields are B = 10 14 G, 2 x 10 14 G, 2.2 x 10 14 G, and 2.3 x 10 14 G from top to 
bottom, respectively. The long dashed lines are the noise amplitude h rrns = vSh (y) for LIGO. 



